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ABSTRACT 



We have studied the galactic cirrus in low surface 
brightness regions using ISOPHOT raster maps at 
90 fim, 150 fim and 180 fim. Observations are used 
to determine dust emission spectra and dust tem- 
peratures. The data extend to longer wavelengths 
than IRAS observation. Compared with DIRBE data 
the resolution and the sensitivity are better and this 
makes it possible to study faint cirrus emission in 
individual fields. 

We will discuss the calibration of the observa- 
tions and present results of a comparison between 
ISOPHOT and DIRBE surface brightness values. 
The correspondence was found to be better than 
~30%. At 90 ^m the ISOPHOT surface brightnesses 
tend to be slightly higher than the DIRBE values 
while at longer wavelength the situation is reversed. 

Surface brightness variations caused by cirrus fluc- 
tuations make it possible to determine the spectrum 
of the dust emission. Cirrus spectra were obtained 
for six fields with surface brightnesses in the range 
l-2MJysr _1 after the Zodiacal light and the contri- 
bution of extragalactic sources have been subtracted. 

Assuming v 2 emissivity law the dust temperatures 
are in the range 18-20 K. Temperature variations can 
be seen even within individual fields. These values 
are higher than values from DIRBE but this can be 
explained by differences in the calibration. 



however, poor spatial resolution and in low surface 
brightness regions the sensitivity of the detectors be- 
comes a limiting factor so that only an average dust 
spectrum over a large area can be determined. 

We use ISOPHOT observations to study the spec- 
trum of galactic cirrus in regions selected for their 
low surface brightness. Because of the better sensi- 
tivity and angular resolution of the ISOPHOT instru- 
ment it is possible to determine the cirrus spectrum 
in these fields based on the cirrus intensity variations. 
In this way constant or slowly varying components, 
e.g. Zodiacal light, can be eliminated. 

The determination of the cirrus spectrum is impor- 
tant also for extragalactic studies where the galac- 
tic foreground emission must be separated. In future 
work the same fields will be used to estimate the FIR 
component of the cosmic background radiation. 

We start by presenting the observations used in this 
study. After that calibration differences between 
ISOPHOT and DIRBE observations are discussed. 
Calibration issues are obviously essential for the de- 
termination of the spectra and the dust tempera- 
tures. Finally, the method used to derive the cirrus 
spectra and the derived dust temperatures are pre- 
sented. A v 2 dust emissivity law will be assumed 
throughout this paper. 



2. OBSERVATIONS 
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1. INTRODUCTION 



At wavelengths longer than ~ 100 fim the far-infrared 
emission of the Galaxy is dominated by thermal 
emission from large dust particles that are in ther- 
mal equilibrium with the interstellar radiation field. 
Dust temperatures determined from DIRBE mea- 
surements are typically 1 6—19 K (e.g. |Lagachc ct al.| 



1998: Dwek et al. 1997) in rough agreement with 



model predictions (e.g. Draine & Anderson 198S) 



Most previous FIR dust temperature determinations 
have been based on COBE observations. These have, 



The cirrus spectra were determined in six fields that 
had been observed at 90 ™. 150 fim and 180 u m with 
the ISOPHOT ins trument ([Lemke ct al. 1996[ ) aboard 
the ISO satellite (Kessler et al. I'JUE). The observa- 
tion consist of PH122 raster maps. Detailed descrip- 
tion of the observations is given in Table n Data 
reduction was done using the IS OPHOT Interactive 
Analysis software package (PIA) ( pabriel et al. 1997] ) 
version 8.0. In the following the analysis is based on 
measurements reduced to the AAP level. 

Flatfielding was done outside PIA using custom rou- 
tines that also removed slow drifts of the detector 
pixels relative to each other. Good flat fielding is not 
critical for the determination of the spectra but it 
will reduce the noise in the subsequent steps. 
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Table 1. Description of the ISOPHOT maps used to determine the cirrus spectra. The columns are: (1) name of 
the field, (2)-(3) coordinates of the centre of each field, (4)-(5) galactic coordinates, (6) map area, (7) ISOPHOT 
filter, (8) number of raster positions, (9) step between adjacent raster positions and (10) integration time. The 
distance of adjacent scans was the same as the raster step used along the scan line. 



Field 


Map Centre 






Area 


Filter 


Rasters 


Step 


tint 




RA(2000.0) 


DEC(2000.0) 


I 


b 


(square degrees) 






(arcsec) 


(") 


(1) 


(2) 


(3) 


(i) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


VCN 


15 15 21.7 


+56 28 58 


91.76 


51.42 


0.030 


C.90 


10x4 


90 


16 














C.135 


10x4 


90 


16 














C.180 


10x4 


90 


46 


VCS 


15 15 53.1 


+56 19 30 


91.27 


51.40 


0.023 


C.90 


21x2 


90 


16 














C.135 


21x2 


90 


16 














C.180 


21x2 


90 


46 


NGPN 


13 43 53.0 


+40 11 35 


86.82 


73.61 


0.27 


C.90 


32x4 


180 


23 














C.135 


32x4 


180 


27 














C.180 


32x4 


180 


27 




13 42 32.0 


+40 29 06 


87.88 


73.26 


0.53 


C.180 


15x15 


180 


32 


NGPS 


13 49 43.7 


+39 07 30 


81.49 


73.30 


0.27 


C.90 


32x4 


180 


23 














C.135 


32x4 


180 


27 














C.180 


32x4 


180 


27 


EBL22 


02 26 34.5 


-25 53 43 


215.78 


-69.19 


0.19 


C.90 


32x3 


180 


23 














C.135 


32x3 


180 


27 














C.180 


32x3 


180 


27 


EBL26 


01 18 14.5 


01 56 40 


135.89 


-60.66 


0.27 


C.90 


32x4 


180 


23 














C.135 


32x4 


180 


23 














C.180 


32x4 


180 


23 



COMPARISON WITH DIRBE SURFACE 
BRIGHTNESS VALUES 



The results presented in this paper are based on the 
calibration performed with Fine Calibration Source 
(FCS) measurements. The absolute calibration of 
ISOPHOT observations has been estimated to be 
be tter than 30 % QKlaas et al. 1998|). [Klaas fc 



Richards (1999) recently reported lor both GlUU and 



C200 accuracies better than 20% relative to DIRBE. 
These results apply to fields with surface brightness 
above ~5MJysr _1 while several of our C200 maps 
are below 4MJysr~ 1 . It is therefore useful to com- 
pare our observations with DIRBE data especially as 
the calibration accuracy is essential for the determi- 
nation of the cirrus spectra. 

The DIRBE and ISOPHOT calibration was studied 
in the case of the fields NGPS, NGPN, EBL22, and 
EBL26 (see Table |l|). Three methods were used to 
perform the comparison: 

1. Comparison of absolute surface brightness lev- 
els. The ISOPHOT maps were compared with 
DIRBE ZSMA (Zodi- Subtracted Mission Aver- 
age) data to which Zodiacal lig ht, was added 



accord ing to the model given by Leinert et al 
(1998)|) . Using the DIRBE weekly maps ob- 
served with the same solar elongation as the 
ISO data did not change the results significantly. 
Because of the lower noise we chose to use the 
ZSMA data. Both DIRBE ZSMA and ISO were 
colour corrected assuming a v 2 B v spectrum with 
Tdust=18K. The zodiacal light estimates were 
colour corrected as suming black body sp e ctrum 



with 7^,^=270 i K ( [Abraham et al. 199 
ham et al. 1999). Corresponding curves 



Abra- 
tted to 

the DIRBE ZSMA 100 ^m, 140 /an and 240 ^m 



data and the zodiacal light estimates were used 
to derive the DIRBE surface brightness esti- 
mates interpolated to the wavelength of the ISO 
observations. The average ISO flux density was 
compared with the derived DIRBE values calcu- 
lated as the average over ISO map and weighted 
by the DIRBE beam. The result is the surface 
brightness ratio S'(DIRBE)/S , (ISOPHOT). 

2. Direct comparison of the surface brightness vari- 
ations in ISOPHOT and DIRBE ZSMA maps. 
All data were first colour corrected for a spec- 
trum v 2 B(y) with T dust =18. OK. For each 
DIRBE pixel the data were interpolated to 
the wavelength of the ISOPHOT observations 
and the corresponding ISOPHOT surface bright- 
ness estimate was calculated as a weighted 
average over the DIRBE beam. Linear fit 
was done to these points to derive the slope 
between the surface brightness values k = 
A5'(DIRBE) / AS*(ISOPHOT) . 

3. Comparison via IRAS. Linear relations were es- 
tablished between the IRAS 100 /mi data and the 
ISO observations and between the IRAS data 
and the DIRBE ZSMA data interpolated to the 
wavelength of the ISO observations using a fit- 
ted curve u 2 B u (T=18K). For each IRAS ISSA 
map pixel inside the field the corresponding ISO 
surface brightness was calculated as an average 
weighted with a Gaussian with FWHM~5 ar- 
cmin. The relation between IRAS and DIRBE 
surface brightnesses was determined with a sim- 
ilar method but over an area with radius ~2 de- 
grees where for each DIRBE pixel the average 
IRAS value was calculated by weighting with 
the DIRBE beam. The slopes of these two lin- 
ear relations AS'(DIRBE)/A5(IRAS(100/im)) 
A5(IRAS(100Atm))/A5(ISOPHOT) give the ra- 
tio between DIRBE and ISO scales, k = 
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AS(DIRBE)/AS(ISOPHOT). 



Method (1) suffers from the small size of the 
ISOPHOT maps which means that for each DIRBE 
pixel a large portion of the flux comes from outside 
the area mapped with ISOPHOT. 

Method (2) is even more affected by the limited map 
sizes and the large size of the DIRBE beam since we 
must determine the relation between brightness vari- 
ations. However, NGPN and NGPS form together 
a 3.1 degrees long strip where the surface brightness 
gradient is along the longer side of the map. Here the 
fact that the ISO map is narrow compared with the 
DIRBE beam should not lead to significant errors. 
On the other hand, in EBL26 the surface bright- 
ness drops quickly outside the bright region in the 
northern part of the map. The surface brightness 
variation is not resolved by DIRBE beam and the 
values A(DIRBE)/A(ISOPHOT) obtained would be 
underestimated. In EBL22 and NGPN the method 
does not work because of the lack of sufficient surface 
brightness variations. 

Methods (2) and (3) are not affected by the presence 
of zodiacal light. The last method uses IRAS ISSA 
maps as an intermediate step and is therefore much 
less affected by the poor resolution of the DIRBE 
data. It allows also the use of DIRBE data over a 
much larger area and thereby reduces the effect of 
the large noise present in 140 /mi and 240 /im DIRBE 
observations. In these low surface brightness regions 
the error estimates of 140 /im and 240 /tm DIRBE 
measurements exceed 50%. On the other hand, the 
method is based on the assumption that the ratio 
of the emission at 100 /im and at other wavelengths 
remains constant. 

The results are given in Table 0. It can be seen 
that DIRBE and ISOPHOT surface brightness val- 
ues agree typically to within 30%. However, at 90 fim 
the ISOPHOT surface brightness tends to be above 
the value predicted from DIRBE observations while 
at longer wavelengths the situation is reversed. The 
extrapolation of the DIRBE values to 90 /im using 
the modified black body fitted to three longer wave- 
length may not be entirely justified and may lead 
to underestimation of the surface brightness. Also, 
in these faintest regions the accuracy of the DIRBE 
140 /im and 240 /tm is poor and a reliable comparison 
with ISOPHOT C200 values is difficult. The modi- 
fied black body curves fitted to DIRBE data are de- 
termined mostly by the 100 /tm DIRBE value which 
has significantly smaller error estimates than the two 
longer wavelengths, and the assumed dust tempera- 
ture, T=18.0K. Increasing the temperature by one 
degree would increase k at 90 /tm by some 5% and 
decrease it at longer wavelengths by some 15%. 

In spite of uncertainties the results suggest that the 
dust temperatures derived from ISOPHOT observa- 
tions are, because of calibration differences, higher 
than those based on DIRBE data. Using PIA ver- 
sions prior to 8.0 the differences between ISOPHOT 
and DIRBE surface brightnesses were typically larger 
by some 20%. 
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Figure 1. Determination of the relative scales of the 
ISOPHOT and DIRBE maps using IRAS 100/dm 
ISSA maps as an intermediate step. The ISOPHOT 
map is the 180 [im map of the field NGPS. The upper 
frame shows linear fit between IRAS 100 (im surface 
brightness values and the DIRBE values interpolated 
to 180 nm. The lower frame shows the correspond- 
ing fit between IRAS 100 [im and ISOPHOT maps. 
The slopes of the fitted lines (shown in figures) can 
be combined to get the value shown in column 5 of 
Table |. 



4. CIRRUS SPECTRA 



As the first step all data were colour corrected in or- 
der to derive monochromatic surface brightness val- 
ues at the reference wavelengths of the filters. Cor- 
rection was done for spectrum v 2 B u (T — 18K) which 
approximates the expected cirrus spectrum. A large 
fraction of the observed surface brightness values is 
due to Zodiacal light. Since we will study only the rel- 
ative surface brightness variations the Zodiacal light 
is eliminated from the analysis and it does not affect 
the required colour correction. 

In the calculations the 180 /tm maps are used as the 
reference. For each 180 /im pixel the corresponding 
surface brightness value at another wavelength is cal- 
culated as a weighted average. The spatial weight- 
ing is done with a Gaussian with FWHMwlOO arcsec 
and the individual measurements are also weighted 
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Table 2. The ratios k between the DIRBE and ISO flux density scales. The columns are: (1) the name of the 
field, (2) wavelength of ISO observations, (3) mean surface brightness of the ISO map, (4) k from the comparison 
of absolute surface brightnesses, (5) k obtained from direct comparison of surface brightness variations and (6) k 
obtained via 100 /im ISSA maps. Formal errors from the least squares fits are shown in parentheses. For the first 
method (i.e. column 4) the numbers given in parentheses give the dispersion in the DIRBE surface brightness 
values interpolated to the wavelength of the ISO observations 



Map 


A 


< S > 


S(DIRBE) /S(ISOPHOT) 


AS(DIRBE) /AS(ISO) 


via IRAS 






(MJv/sr) 








(1) 

V- 1 -/ 


(2) 


(3) 




(5) 


(6) 


EBL22 


90 


5.6 


0.81(0.01) 








150 


3.4 


1.28(0.03) 








180 


2.V 


1.42(0.05) 






EBL26 


00 


14.3 


0.84(0.02) 




0.88(0.10) 




150 


9.1 


1.04(0.05) 




1.09(0.12) 




180 


7.4 


1.09(0.03) 




1.27(0.14) 


NGPS 


00 


5.7 


0.69(0.06) 


0.69(0.10) 


0.70(0.05) 




150 


3.7 


0.94(0.15) 


1.68(0.15) 


1.43(0.11) 




180 


3.1 


1.04(0.02) 


1.59(0.15) 


1.40(0.10) 


NGPN 


90 


5.0 


0.67(0.02) 








150 


3.6 


0.97(0.04) 








180 


2.7 


1.19(0.11) 








180 1 


2.6 


1.21(0.07) 






NGPN & NGPS 


90 


5.2 


0.67(0.02) 


0.63(0.21) 


0.75(0.06) 




150 


4.1 


0.92(0.04) 


0.91(0.34) 


1.56(0.13) 




180 


3.3 


1.08(0.12) 


0.84(0.33) 


1.53(0.11) 



1 the larger 180 fira map 



ac( ording to their error estimates. The calculated 
values are plotted against the 180 /im surface bright- 
ness and a straight line is fitted to the points with 
a least squares algorithm that takes into account the 
error estimates on both axes. The procedure is re- 
peated for all maps and the slopes of the fitted lines 
provide the emission spectrum i.e. for each observed 
wavelength the relative intensity with respect to the 
180 /im emission. 

Figure ^ shows as an example the relations between 
180 /im and 90 /im and 180 /im and 150 /im surface 
brightness values in the field EBL26. The fitted least 
square lines are also shown. Zodiacal light will ap- 
pear merely as an offset and does not affect the slopes 
that determine the cirrus spectrum. Furthermore, zo- 
diacal light is weak at 180 /im and we are looking only 
for variations correlated with 180 /xm emission. 

Galactic cirrus is the dominant source of brightness 
fluctuations and therefore the spectra obtained char- 
acterize the galactic cirrus emission. However, in low 
surface brightness regions there may be a significant 
contribution from faint extragalactic sources. The 
DIRBE ( |Hauser et al. 1998j ) and FIRAS flFixen et al 
19^8|) experiments indicated a FIR cosmic infrared 
background flux of 1 MJy sr _1 between 100 and 
240 /im. In fields like EBL22 this would mean that 
at 180 /im about one third of the surface brightness 
is due to extragalactic sources. Since the cirrus spec- 
tra are determined from the brightness variations we 
must consider the brightness fluctuations caused by 
these sources in relation to the cirrus fluctuations. 

The cirrus power spectrum is proportional to Bo/k 3 , 
where B n is the cirrus brightn ess and k the spatial fre- 
quency ( Gautier et al. 1992 ). In other words, cirrus 
fluctuations decrease rapidly as we mov e to fainter 
cirrus regions and smaller spatial scales. Lagache 



Puget (1999) have reported a detection of the ex- 
tragalactic background fluctuations at 175 /im in the 
Marano 1 field that has cirrus surface brightness com- 
parable to our fields. Even after removal of detected 
sources they found at scales below 10 arcmin fluctu- 
ations in excess of the expected cirrus contribution. 
This was interpreted to be caused by point sources 
below the detection limit. 

The larger dimension of our maps is typically 1.5 de- 
grees and most of the surface brightness variations 
can be expected to be due to cirrus. This is particu- 
larly clear in fields NGPS and EBL26 where there is 
a clear intensity gradient along the mapped strip con- 
sistent with the larger scale cirrus distribution visible 
in IRAS maps. On the other hand, the field EBL22 
is remarkably flat and it is conceivable that in the 
absence of large cirrus structures the extragalactic 
point sources could affect the derived spectrum. The 
small dynamic range means, however, that in this 
case the accuracy of the derived spectrum is not very 
good. The fields VCN and VCS were observed close 
the bright galaxy NGC 5907 and the VCS field ac- 
tually extends over the centre of the galaxy. Pix- 
els that were clearly affected by the emission from 
NGC 5907 were first discarded but it is still possible 
that the derived spectra are affected by emission of 
the galaxy. Because of the small field size (~10 ar- 
cmin) the VCN spectrum could be affected by faint 
extragalactic sources. 

The spectra obtained for the six fields are presented 
in Figure a. The error bars include the error esti- 
mates obtained from the fitting procedure and a fur- 
ther 10% error representing the error in the relative 
calibration of the observations made with different 
filters. Modified Planck curves, v 2 B v {T\ were fitted 
to the data and the dust temperatures obtained are 
shown in the figure. 
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Figure 2. Plots of the 90 fim and the 150 [im sur- 
face brightness values in the EBL26 field as function 
of 180 \ira surface brightness. The fitted least squares 
lines are also shown. The slopes of these lines give the 
ratio between surface brightnesses at different wave- 
lengths i. e. the spectrum of the cirrus emission. 



5. DISCUSSION 



We have determined cirrus dust temperatures in six 
fields with low surface brightness. The derived dust 
temperatures, ~18-20K, are higher than the typi- 
cal temperatures estimated for galactic dust from 
DIRBE and FIRAS experiments assuming the same 
v 2 emissivity law. However, the values are still 
within the r ange found in these studies. Lagache 



et al. (1998)| found an average temperature ot 17.5 K 
while higher temperatures were mostly associated 
with star forming regions. Temperatures exceeding 
20 K were seen only towards the galactic centre and 
some prominent star forming regions although the 
scatter in the temperatures is more than 1 K at all 
galactic latitudes. 

Our fields are all at high galactic latitudes (see Ta- 
ble [j]). Based on COBE measurements at 100-240 /im 
Boulanger et al. (1996) derived average dust temper- 
at ures of T = 17.5 K for regions at | b | >30 degrees. 



Earlier FIRAS estimates (pleach et al. 1995) also sug- 
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Figure 3. Cirrus spectra derived from the six observed 
fields. The error estimates contain the statistical er- 
rors provided by the fitting procedure and a further 
10% error representing the uncertainty of the rela- 
tive calibration between the filters. The solid lines 
show the fitted modified black body curves, v 2 B„{T). 



with |6| >45 degrees. 

In Section bj we found systematic differences be- 
tween ISOPHOT and DIRBE calibrations. Com- 
pared with DIRBE the ISOPHOT surface bright- 
nesses were found to be higher at 90 /im and lower 
at longer wavelengths. According to Table || adopt- 
ing the DIRBE calibration would increase our 90 /im 
data by up to 30% relative to the longer wavelengths. 
This would reduce the estimated dust temperature 
by ^1.5 K. Therefore the temperature difference be- 
tween DIRBE and our results can be explained en- 
tirely by the difference in calibration. 

Our temperature determinations are based on data at 
three wavelengths: 90 /im, 150 /im and 240 /im. The 
emission at wavelengths A >100/im can be explained 
by classical grains in equilibrium with the interstellar 
radiation field. Smaller grains that are transiently 
heated to higher temperatures are responsible for 
the emission at shorter wavelengths. Based on both 
observations and theoretical calculations the contri- 
bution from transietly heated particles is significant 
already at 60 /im a nd may extend even to slightly 



gested similar or lower values. H owever, based on the 



longer wavelengths (Sodroski et al. 1994 ; Dwek et al 
1997| ; [Desert et al. P)90| ). The ISOPHOT 90 /mi filter' 



DIRBE 140 /im to 240 /im ratios jDwck et al. (1997)| ) 
reported colour temperatures of T ~19 K tor regions 



is wide enough to pick up some radiation from wave- 
lcnths down to ^60 fim and our observations may be 
affected to a small degree by this other grain popula- 
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tion. Our temperature estimates do not correspond 
exactly to the equilibrium temperature of large dust 
particles but the bias towards higher temperatures is 
expected to be small. 

In the low surface brightness regions the accuracy 
of DIRBE measurements is not very good and one 
can obtain only the average spectrum of a very large 
area. Because of the better resolution and sensitivity 
of ISOPHOT we have been able to study the cirrus 
spectrum in separate, small regions. Adopting the 
ISOPHOT calibration the accuracy of the tempera- 
ture estimates is ~1K. For example, the 2K differ- 
ence between the adjacent fields NGPS and NGPN 
is likely to be true. 
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